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This chapter presents an introduction to graph colouring algorithms^ The fo¬ 
cus is on vertex-colouring algorithms that work for general classes of graphs 
with worst-case performance guarantees in a sequential model of computa¬ 
tion. The presentation aims to demonstrate the breadth of available tech¬ 
niques and is organized by algorithmic paradigm. 


1. Introduction 

A straightforward algorithm for finding a vertex-colouring of a graph is to search 
systematically among all mappings from the set of vertices to the set of colours, a 
technique often called exhaustive or brute force: 

Algorithm X {Exhaustive search) Given an integer q > \ and a graph G with 
vertex set V, this algorithm finds a vertex-colouring using q colours if one exists. 

XI [Main loop] For each mapping /: V —> {1,2,..., < 7 ), do Step X2. 

X2 [Check /] If every edge vw satisfies /(v) t f{w), ferminafe wifh / as fhe 
resulf. ■ 

This algorifhm has few redeeming qualifies, ofher fhan ifs being correcf. We 
consider if here because if serves as an opporfunify fo make explicif fhe framework 
in which we presenf more inferesfing algorifhms. 

’ Appears as Thore Husfeldt, Graph colouring algorithms. Chapter XIII of Topics in Chromatic 
Graph Theory, L. W. Beineke and Robin J. Wilson (eds.), Encyclopedia of Mathematics and its 
Applications, Cambridge University Press, ISBN 978-I-I07-03350-4, 2015, pp. 277-303. 
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Model of computation 


If G has n vertices and m edges, then the number of operations used by Algorithm 
X can be asymptotically bounded by 0{q'^{n + m)), which we call the running time 
of the algorithm. 

To make such a claim, we tacitly assume a computational model that includes 
primitive operations, such as iterating over all mappings from one finite set A to 
another finite set B in time (Step XI), or iterating over all edges in time 

0{n + m) (Step X2). For instance, we assume that the input graph is represented by 
an array of sequences indexed by vertices; the sequence stored at vertex v contains 
the neighouring vertices N{v), see Fig. [T] This representation allows us to iterate 




Fig. 1: A graph and its representation as an array of sequences 

over the neighbours of a vertex in time 0(deg v). (An alternative representation, 
such as an incidence or adjacency matrix, would not allow this.) Note that detecting 
whether two graphs are isomorphic is not a primitive operation. The convention of 
expressing computational resources using asymptotic notation is consistent with 
our somewhat cavalier attitude towards the details of our computational model. 
Our assumptions are consistent with the behaviour of a modem computer in a high- 
level programming language. Nevertheless, we will explain our algorithms in plain 
English. 


Worst-case asymptotic analysis 

Note that we could have fixed fhe colouring of a specific verfex v as /(v) = 0, 
reducing Algorifhm X’s running time fo 0{q''~^{n + m)). A momenf’s fhoughf 
shows fhaf fhis reasoning can fhen be exfended fo cliques of size r > 1: search 
fhrough all induced subgraphs until a clique of size r is found, arbifrarily map 
fhese vertices fo {1,2,..., r) and fhen lef Algorifhm X colour fhe remaining verfices. 
This reduces fhe running lime fo 0{q'^~‘^^^^n‘^^^\n + m)), where a»(G) is fhe clique 
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size. This may be quite useful for some graphs. Another observation is that in the 
best case, the running time is 0{n + m). However, we will normally not pursue 
this kind of argument. Instead, we are maximally pessimistic about the input and 
the algorithm’s underspecified choices. In other words, we understand running 
times as worst-case performance guarantees, rather than ‘typical’ running times or 
average running times over some distribution. 

Sometimes we may even say that Algorithm X requires time q” poly(n), where 
we leave the polynomial factor unspecified in order fo signal fhe perfunctory affen- 
fion we extend to fhese issues. 

Overview and notation 

Straightforward variants of Algorithm X can be used to solve some other graph 
colouring problems. For instance, to find a list-colouring, we restrict the range of 
values for each /(v) to a given list; to find an edge-colouring, we iterate over all 
mappings /: E {1,1,... ,q}. 

Another modification is to count the number of colourings instead of finding 
just one. These extensions provide baseline algorithms for list-colouring, edge¬ 
colouring, the chromatic polynomial, the chromatic index, and so forth. However, 
for purposes of exposition, we present algorithms in their least general form, em¬ 
phasizing the algorithmic idea rather than its (sometimes quite pedestrian) general¬ 
izations. The algorithms are organized by algorithmic technique rather than prob¬ 
lem type, graph class, optimality criterion, or computational complexity. These 
sections are largely independent and can be read in any order, except perhaps for 
Algorithm G in Section|2l The final section takes a step back and relates the various 
colouring problems to each other. 

2. Greedy colouring 

The following algorithm, sometimes called the greedy or sequential algorithm, con¬ 
siders the vertices one by one and uses the first available colour. 

Algorithm G {Greedy vertex-colouring) Given a graph G with maximum degree A 
and an ordering vi, V 2 ,..., v„ of its vertices, this algorithm finds a vertex-colouring 
with max, \{ j < i: vjvt e£')|-i-l <A-i-l colours. 

G1 [Initialize] Set t = 0. 

G2 [Next vertex] Increment i.lli - n -v \, terminate with / as the result. 

G3 [Find the colours N(v,)] Compute the set C = \Jj<if{vj) of colours already 
assigned to the neighbours of v,-. 
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G4 [Assign the smallest available colour to v,] For increasing c = 1,2,, check 
whether c e C. If not, set /(v,) = c and return to Step G2. ■ 

For the number of colours, it is clear that in Step G4, the value of c is at most 
|C|, which is bounded by the number of neighbours of v,- among vi, V2, ..., v,_i. In 
particular. Algorithm G establishes that;^(G) < A(G) + 1. 

For the running time, note that both Steps G3 and G4 take at most 0(1 + deg v,) 
operations. Summing over all i, the total time spent in Steps G3 and G4 is asymp¬ 
totically bounded by n -i- (deg vi + deg V2 -i- 1 - deg v„) = n + 2m. Thus, Algorithm 

G takes time 0{n + m). 

Optimal ordering The size of the colouring computed by Algorithm G depends 
heavily on the vertex ordering. Its worst-case behaviour is poor. For instance, it 
spends colours on the 2-colourable crown graph shown in Fig.|2j 


Vl V3 V5 V„_i 



V2 V4 V6 V„ 


Fig. 2: The crown graph 

On the other hand, for every graph there exists an ordering for which Algo¬ 
rithm G uses an optimal number of colours; indeed, any ordering that satisfies 
fiVi) < f{Vi+0 for an optimal colouring / has this property. Since there are n\ 
different orderings, this observation is algorithmically quite useless. An ordering 
is perfect for a graph if, for every induced subgraph. Algorithm G results in an op¬ 
timal colouring; triangulated graphs and comparability graphs always admit such 
an ordering, as shown by Chvatal ifTTl . 

Randomness 

Algorithm G performs quite well on random graphs, whatever the vertex order¬ 
ing. For almost all n-vertex graphs, it uses n/(log n-3 log log n) colours, which is 
roughly twice the optimum value (see liTSll '). 

This suggests the following randomized algorithm. For a graph G, choose a 
vertex ordering at random and then execute Algorithm G. For many problems, it is 
a sound algorithmic design strategy to trade good average-case behaviour for good 
(expected) worst-case behaviour in this way. However, for Algorithm G the result 
is quite poor: for every e > 0 there exist graphs with chromatic number n‘^ for 
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which the randomized algorithm uses Q{n/ log n) colours with high probability, as 
shown by Kucera i26l . 


Other orderings 

In the largest-first vertex-degree ordering introduced by Welsh and Powell |[38l . 
the vertices are ordered such that degvi > degV 2 > • • • > degv„. This establishes 
the bound;t^(G) < 1 -i- max,- min{deg v,-, i - 1), which is sometimes better than 1 -i- A, 
such as in Fig. |3] 



Fig. 3: 

Closely related in spirit is Matula’s smallest-last ordering Il32ll . given as fol¬ 
lows: choose as the last vertex v„ a vertex of minimum degree in G, and proceed 
recursively with G - v„, see Fig. ID With this ordering, the size of the resulting 
colouring is be bounded by the Szekeres-Wilf bound 1361, 

XiG) < dgn(G) + 1, 

where the degeneracy dgn(G) is the maximum over all subgraphs H of G of the 
minimum degree S{H). This ordering optimally colours crown graphs and many 
other classes of graphs, and uses six colours on any planar graph. 

Vg Vi V4 V3 V2 Vj vg V4 V2 Vi V3 V5 

largest-first: . smallest-last: . 


Fig. 4: 

Other orderings are dynamic in the sense that the ordering is determined during 
the execution of the algorithm, rather than in advance. For example, Brelaz ||6l sug¬ 
gests choosing the next vertex from among those adjacent to the largest number of 
different colours. Many other orderings have been investigated (see the surveys of 
Kosowski and Manuszewski Ii25l and Maffray lOTl l. Many of them perform quite 
well on instances that one may encounter ‘in practice’, but attempts at formalizing 
what this means are quixotic. 
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2-colourable graphs 

Of particular interest are those vertex orderings in which every vertex v, is adjacent 
to some vertex vj with j < i. Such orderings can be computed in time 0{m + n) 
using basic graph-traversal algorithms. This algorithm is sufficiently important to 
be made explicit. 

Algorithm B {Bipartition) Given a connected graph G, this algorithm finds a 
2-colouring if one exists. Otherwise, it outputs an odd cycle. 

B1 [Initialize] Let /(vi) = 1 and let Q (the ‘queue’) be an empty sequence. For 
each neighbour w of vi, set p{w) = vi (the ‘parent’ of w) and add w to Q. 

B2 [Next vertex] If Q is empty, go to Step B3. Otherwise, remove the first vertex 
V from Q and set /(v) to the colour not already assigned to p{v). For each 
neighbour w of v, if w is not yet coloured and does not belong to Q, then set 
p(w) - V and add w to the end of Q. Repeat Step B2. 

B3 [Verify 2-colouring] Iterate over all edges to verify that /(v) f{w) for every 
edge vw. If so, terminate with / as the result. 

B4 [Construct odd cycle] Let v>v be an edge with /(v) = /(w) and let u be the 
nearest common ancestor of v and w in the tree defined by p. Output the path 
w, p{w), p{p{w)),.. .,u, followed by the reversal of the path v, p{v), p{p{v)), 
... ,u, followed by the edge v>v. ■ 


Fig. |5] shows an execution of Algorithm B finding a 2-colouring. 



Fig. 5: Execution of Algorithm B 


Algorithm B is an example of a ‘certifying’ algorithm: an algorithm that pro¬ 
duces a witness to certify its correctness, in this case an odd cycle if the graph is 
not 2-colourable. To see that the cycle constructed in Step B4 has odd length, note 
that on the two paths w, p(w), p{p{w)), ...,u and v, p{v), p{p{v)),... ,u, each vertex 
has a different colour from its predecessor. Since the respective endpoints of both 
paths have the same colour, they must contain the same number of edges modulo 
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2. In particular, their total length is even. With the additional edge vw, the length 
of the resulting cycle is odd. 

The order in which the vertices are considered by Algorithm B depends on the 
first-in first-out behaviour of the queue Q. The resulting ordering is called breadth- 
first. An important variant uses a last-in first-out ‘stack’ instead of a queue; the 
resulting ordering is called depth-first. Fig. 0 shows the resulting behaviour on the 
graph from Fig. [51 



Fig. 6: 



Execution of Algorithm B using depth-first search 


Algorithm B works also for the list-colouring problem, provided that for each 
vertex v, the available list of colours L(v) has size at most 2. This observation leads 
to a simple, randomized, exponential-time algorithm for 3-colouring due to Beigel 
and Eppstein [Tj 

Algorithm P {Palette restriction) Given a graph, this algorithm finds a 3-colouring 
if one exists. 

PI [Eorbid one colour at each vertex] Eor each vertex v, select a list L(v) of colours 
available at v uniformly and independently at random from the three lists {1,2), 
{2,3), and {1,3). 

P2 [Attempt 2-colouring] Try to solve the list-colouring instance given by L using 
Algorithm B, setting f{v\) - minL(vi) in Step BE If successful, terminate 
with the resulting colouring. Otherwise, return to Step PI. ■ 

To analyse the running time, consider a 3-colouring /. Eor each vertex v, colour 
/(v) belongs to L(v) with probablity |. Thus, with probability at least (|)”, the list 
colouring instance constructed in step PI has a solution. It follows that the expected 
number of repetitions is (|)", each of which takes polynomial time. 

Wigderson’s algorithm 

Algorithms B and G appear together in Wigderson’s algorithm liTOl : 







CHAPTER 13. GRAPH COLOURING ALGORITHMS 


Algorithm W {Wigderson’s algorithm) Given a 3-chromatic graph G, this algo¬ 
rithm finds a vertex-colouring with 0{s/n) colours. 

W1 [Initialize] Letc = 1. 

W2 [A(G) — rV^l] Consider a vertex v in G with degv > if such ver¬ 
tex exists, go to Step W3. Use Algorithm B to 2-colour the neighbourhood 
G[A(v)] with colours c and c -i- 1. Remove N{v) from G and increase c by 
;y(G[A(v)]). Repeat Step W2. 

W3 [A(G) < rV^l] Use Algorithm G to colour the remaining vertices with the 

colours c, c-I-1,..., c-I-[V^l- B 

Fig. |7] shows an execution of Algorithm W finding a 5-colouring of the 16- 
vertex instance from Fig. [T] 



Fig. 7: Execution of Algorithm W 

The running time is clearly bounded by 0{n + m). To analyse the number of 
colours, we first need to verify Step W2. Since G is 3-colourable, so is the subgraph 
induced by N{v) U (v). Now, if G[A(v)] requires 3 colours, then G[N{v) U {v}] 
requires 4, so G[A(v)] is 2-colourable and therefore Step W2 is correct. Note that 
Step W2 can be run at most 0(^Jn) times, each using at most two colours. Step W3 
expends another [ sjn) colours according to Algorithm G. 

Algorithm W naturally extends to graphs with ;y(G) > 3. In this case. Step 
W2 calls Algorithm W recursively to colour (;y(G) - l)-colourable neighbourhoods. 
The resulting algorithm uses G(n^“F(i-^(G))) colours. 


3. Recursion 

Recursion is a fundamental algorithmic design technique. The idea is to reduce a 
problem to one or more simpler instances of the same problem. 
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Contraction 

The oldest recursive construction for graph colouring expresses the chromatic poly¬ 
nomial P{G, q) and the chromatic number ;y(G) in terms of edge-contractions: For 
non-adjacent vertices v, w and integer q - 9,1, ■ ■ ■ ,n, 

P{G, q) = P{G -F vw, q) + P(G/v>v, q) , 

X{G) = minixiG + vw),x{Glvw)], 

see Chapter 3, Section 2.1. These ‘addition-contraction’ recurrences immediately 
imply a recursive algorithm. For instance, 

= PiK4, q) + P(^3, -?) - q{q - 1)(-? - 2){{q - 3){q - 4) + 1) . 

Note that the graphs at the end of the recursion are complete. 

For sparse graphs, it is more useful to express the same idea as a ‘deletion- 
contraction’ recurrence, which deletes and contracts edges until the graph is empty: 

P(G, q) = P{G/e, q) - P{G - e, q) (e e E). 

Many other graph problems beside colouring can be expressed by a deletion- 
contraction recurrence. The most general graph invariant that can be defined in this 
fashion is the Tutte polynomial (see iS and ifT^ for its algorithmic aspects). 

The algorithm implied by these recursions is sometimes called Zykov’s algo¬ 
rithm ll42l . Here is the deletion-contraction version. 

Algorithm C {Contraction) Given a graph G, this algorithm returns the sequence 
of coefficients (oq, ai,..an) of the chromatic polynomial P{G, q) = X"=o 

Cl [Base] If G has no edges then return the coefficients (0,0,..., 0,1), correspond¬ 
ing to the polynomial P{G, q) = q^. 

C2 [Recursion] Pick an edge e and construct the graphs G' = Gje and G" - 
G - e. Call Algorithm C recursively to compute P{G',q) and P{G",q) as 
sequences of coefficients (a^, rz',..., and (a”, a"a"). Return {a'^ - 
a'^',a^ - a'f ... ,a'n- a”), corresponding to the polynomial P{Gle, q) - P{G - 
e,q). ■ 

To analyse the running time, let T{r) be the number of executions of Step Cl for 
graphs with n vertices and m edges, where r - n + m. The two graphs constructed 
in Step C2 have size n - I + m - 1 = r - 2 and n + m - 1 = r - I, respectively. 
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so T satisfies T{r) = T{r - \) + T{r - 2). This is a well-known recurrence with 
solution T{r) - where (p = ^{\ + yJS) is the golden ratio. Thus, Algorithm 

C requires poly(n) = 0(1.619"’'''") time. A similar analysis for the algorithm 
implied by the deletion-addition recursion gives poly(n), where rn = - m 

is the number of edges in the complement of G. 

These worst-case bounds are often very pessimistic. They do not take into 
account that recurrences can be stopped as soon as the graph is a tree (or some 
other easily recognized graph whose chromatic polynomial is known as a closed 
formula), or that P factorizes over connected components. Moreover, we can use 
graph isomorphism heuristics and tabulation to avoid some unnecessary recompu¬ 
tation of isomorphic subproblems (see Ull). Thus, Algorithm C is a more useful 
algorithm than its exponential running time may indicate. 

Vertex partitions and dynamic programmin g 

We turn to a different recurrence, which expresses xiG) in terms of induced sub¬ 
graphs of G. By taking S to be a colour class of an optimal colouring of G, we 
observe that every graph has an independent set of vertices S for which ;y(G) = 
1 +xiG - S). Thus, we have 

XiG) = l+mmxiG-S), (1) 

where the minimum is taken over all non-empty independent sets 5 in G. 

The recursive algorithm implied by ([T]) is too slow to be of interest. We expedite 
it using the fundamental algorithmic idea of dynamic programming. The central 
observation is that the subproblems ;y(G-5) for various vertex-subsets S appearing 
in ([T]) are computed over and over again. It thus makes sense to store these 2" 
values in a table when they are first computed. Subsequent evaluations can then be 
handled by consulting the table. 

We express the resulting algorithm in a bottom-up fashion: 

Algorithm D {Dynamic programming) Given a graph G, this algorithm computes 
a table T with T{W) - x{G[W]), for each W QV. 

D1 [Initialize] Construct a table with (initially undefined) entries T{W) for each 
W c V. Set r(0) - 0. 

D2 [Main loop] List all vertex-subsets Wi,W 2 ,... ,W 2 " £ V in non-decreasing 
order of their size. Do Step D3 for W = W 2 , W 3 ,..., W 2 ", then terminate. 

D3 [Determine r(W)] Set T{W) = 1 - 1 - min T{W \S), where the minimum is taken 
over all non-empty independent sets S in G[W]. ■ 
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The ordering of subsets in the main loop D2 ensures that each set is handled 
before any of its supersets. In particular, all values T{W\S) needed in Step D3 will 
have been previously computed, so the algorithm is well defined. The minimization 
in Step D3 is implemented by iterating over all 2^^^ subsets of W. Thus, the total 
running time of Algorithm D is within a polynomial factor of 



( 2 ) 


This rather straightforward application of dynamic programming already pro¬ 
vides the non-trivial insight that the chromatic number can be computed in time 
exponential in the number of vertices, rather than depending exponentially on m, 
X{G), or a superlinear function of n. 


Maximal independent sets 


To pursue this idea a little further we notice that S in ([I]) can be assumed to be a 
maximal independent set - that is, not a proper subset of another independent set. 
To see this, let / be an optimal colouring and consider the colour class S - /“^(l). 
If S is not maximal, then repeatedly pick a vertex v that is not adjacent to S, and 
set/(v) ^ 1. 

By considering the disjoint union of triangles, we see that there exist k- 
vertex graphs with 3^^^ maximal independent sets. It is known that this is also an 
upper bound, and that the maximal independent sets can be enumerated within a 
polynomial factor of that bound (see fT*], |[341 and OTl l. We therefore have the 
following result: 

Theorem 3.1 The maximal independent sets of a graph on k vertices can be listed 
in time 0(3^^^) and polynomial space. 

We can apply this idea to Algorithm D. The minimization in Step D3 now takes 
the following form: 

D3' [Determine r(W)] Set T{W) = 1 -i- min r(W\S), where the minimum is taken 
over all maximal independent sets S in G[VT]. 

Using Theorem 13. II with k = |W| for the minimization in Step D3', the total 
running time of Algorithm D comes within a polynomial factor of 



For many years, this was the fastest known algorithm for the chromatic number. 
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3-colouring 

Of particular interest is the 3-colouring case. Here, it makes more sense to let 
the outer loop iterate over all maximal independent sets and check whether the 
complement is bipartite. 

Algorithm L {Lawler’s algorithm) Given a graph G, this algorithm finds a 3- 
colouring if one exists. 

LI [Main loop] For each maximal independent set S of G, do Step L2. 

L2 [Try f{S ) = 3] Use Algorithm B to find a colouring / : V\S ^ {1,2} of G-S 
if one exists. In that case, extend / to all of V by setting f(v) - 3 for each 
V e S, and terminate with / as the result. ■ 

The running time of Algorithm L is dominated by the number of executions 
of L2, which, according to Theorem 13.11 is 3”^^. Thus, Algorithm L decides 3- 
colourability in time 3"^^ poly(n) = 0(1.442”) and polynomial space. 

The use of maximal independent sets goes back to Christofides ifTOll . while 
Algorithms D and L are due to Lawler |[28l . A series of improvements to these 
ideas have further reduced these running times. At the time of writing, the best- 
known time bound for 3-colouring is 0(1.329”) by Beigel and Eppstein |T|. 

4. Subgraph expansion 

The Whitney expansion |[39ll of the chromatic polynomial is 

F(G,^) = 2](-1)I^I/(^>; 

ACE 

see Chapter 3, Section 2 for a proof. It expresses the chromatic polynomial as 
an alternating sum of terms, each of which depends on the number of connected 
components k{A) of the edge-subset A Q E. Determining k{A) is a well-studied 
algorithmic graph problem, which can be solved in time 0{n + m) (for example, by 
depth-first search). Thus, the Whitney expansion can be evaluated in time 0{2'"{n + 
m)). 


A more recent expression (see Ej) provides an expansion over induced sub¬ 
graphs: 
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Theorem 4.1 For W C V, let g(VK) be the number of non-empty independent sets 
in G\W]. Then G can be q-coloured if and only if 

WCV 

Proof For each W c V, the term (g(lF))^ counts the number of ways of selecting 
q non-empty independent sets 5'i,5'2,.. .,Sq, where Si c W. For U Q V, let h{U) 
be the number of ways of selecting q non-empty independent sets whose union is 
U. Then = T.ucwKU), so 

WCV WCV ucw 

UCV WDU 

For the last step, note that the inner sum (over IF, with U CW Q V) vanishes except 
when U = V, because there are as many odd-sized as even-sized sets sandwiched 
between different sets, by the principle of inclusion-exclusion. 

If h{V) is non-zero, then there exist independent sets S i,S 2 , ■ ■ ■ ,Sq whose 
union is V. These sets correspond to a colouring: associate a colour with the 
vertices in each set, breaking ties arbitrarily. ■ 

For each IF c F, we can compute the value g(lF) in time 0{2^^^m) by con¬ 
structing each non-empty subset of IF and testing it for independence. Thus, the 
total running time for evaluating Q is within a polynomial factor of 3", just as in 
the analysis ^ for Algorithm D; however, the space requirement here is only poly¬ 
nomial. We can further reduce the running time to 0(2.247”) by using dedicated 
algorithms for evaluating g(lF) from the literature (see Q). 

If exponential space is available, we can do even better. To that end, we first 
introduce a recurrence for g. 

Theorem 4.2 Let IF c F. VFe have g{<Z>) = 0, and, for every v € IF, 

g{W) - g(lF \ {v}) + g(lF \ N[v]) + 1. (4) 

Proof Fix V € IF. The non-empty independent sets 5 c IF can be partitioned into 
two classes with v ^ S and v € S. In the first case, S is a non-empty independent 
set with S c IF \ (v) and thus accounted for by the first term of dUl. Consider the 
second case. Since S contains v and is independent, it contains no vertex from N{v). 
Thus, 5 is a non-empty independent set with (v) c S c IF \ N{v). The number 
of such sets is the same as the number of (not necessarily non-empty) independent 
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sets S' with S' C W \ A^[v], because of the bijective mapping S S' where 
S' = S \ {v). By induction, the number of such sets is giyV \ A^[v]) + 1, where the 
‘+1 ’ term accounts for the empty set. ■ 

This leads to the following algorithm, due to Bjorklund et al. [J] : 


Algorithm I (Inclusion-exclusion) Given a graph G and an integer q > \, this 
algorithm determines whether G can be q-coloured. 

11 [Tabulate g] Set g(<Z) = 0. For each non-empty subset IT c F in inclusion 

order, pick v € IT and set g(lT) = g{W \ {v}) -i- g(W \ A[v]) -i- 1. 

12 [Evaluate ©J If > 0 output ‘yes’, otherwise ‘no’. ■ 


Both Steps II and 12 take time 2”poly(n), and the algorithm requires a table 
with 2” entries. Fig.|8]shows the computations of Algorithm I on a small graph for 
q - 2 and q - 3, with aq(W) = (-l)l^'''^'(g(lT))^. The sum of the entries in column 
a 2 is 0, so there is no 2-colouring. The sum of the entries in column aj, is 18, so a 
3-colouring exists. 
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Fig. 8: Execution of Algorithm I 


With slight modifications. Algorithm I can be made to work for other colouring 
problems such as the chromatic polynomial and list-colouring, also in time and 
space 2”poly(n) (see [31); currently, this is the fastest known algorithm for these 
problems. For the chromatic polynomial, the space requirement can be reduced to 
0(1.292”), while maintaining the 2”poly(?i) running time (see [Ti]). 

5. Local augmentation 

Sometimes, a non-optimal colouring can be improved by a local change that re¬ 
colours some vertices. This general idea is the basis of many local search heuristics 
and also several central theorems. 
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Kempe changes 


An important example, for edge-colouring, establishes Vizing’s theorem, A(G) < 
x'{G) < A(G) -I- 1. Chapter 5 gives a modern and more general presentation of 
the underlying idea, and our focus in the present chapter is to make the algorithm 
explicit. 

A colour free at v if it does not appear on an edge at v. (We consider an 
edge-colouring with A(G) -i-1 colours, so every vertex has at least one free colour.) 
A (Vizing)/a?i around v is a maximal set of edges vwq, vwi,..., vwy, where vwq is 
not yet coloured and the other edges are coloured as follows. For j = 0,1,..., r, 
no colour is free at both v and wj. For j = 1,2,..., r, the jth fan edge vwj has 
colour j and the colours appearing around wj include 1,2 ,..., y but not y -i- 1 (see 
Fig. |9ta)). Such a fan allows a recolouring by moving colours as follows: remove 
the colour from vwy and set /(vwq) = l,/(vwi) = 2,... ,/(vwy-i) = y. This is 
called downshifting from j (see Fig.|9tb)). 



W2 Wj 


Wo 




Wr 


V r + 1 


V .A 




(c) 


Fig. 9: (a) A fan (b) Downshifting from 3 (c) Step V7: colour y is free at Wr+\ 


Algorithm V {Vizing’s algorithm) Given a graph G, this algorithm finds an edge 

colouring with at most A{G) -l- 1 colours in time 0{nm). 

VI [Initialize] Order the edges arbitrarily e\,e 2 , ■ ■ ■ ,e,n. Let i = 0. 

V2 [Extend colouring to next edge] Increment i. If i = ni + I then terminate. 
Otherwise, let vw = e/. 

V3 [Easy case] If a colour c is free at both v and w, then set f(vw) = c and return 
to Step V2. 

V4 [Eind wq and wi] Eet wq = w. Pick a free colour at wq and call it 1. Eet vwi be 
the edge incident with v coloured 1. (Such an edge exists because 1 is not also 
free at v.) 
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V5 [Find W 2 \ Pick a free colour at wi and call it 2. If 2 is also free at v, then set 
/(vivo) - Ij f{vw\) = 2, and return to Step V2. Otherwise, let vw 2 be the edge 
incident with v coloured 2. Set r = 2. 

V6 [Extend fan to Wr+i] Pick a free colour at Wr and call it r + 1. If r + 1 is also 
free at v then downshift from r, recolour f(vwr) - Cr+i and return to Step V2. 
Otherwise, let vw^+i be the edge incident with v coloured r + 1. If each colour 
1,2,..., r appears around Wr+i, then increment r and repeat Step V6. 

V7 [Build a [0,/-path from wj or from Wr+i] Let j € [1,2,..., r) be a free colour 
at Wr+i and let 0 be a colour free at v and different from j. Construct two 
maximal [0, /-coloured paths Pj and P^+i from wj and iv^+i, respectively, by 
following edges of alternating colours 0, j, 0, / ... (see Fig. He)). (The paths 
cannot both end in v.) Let k = j or r + I so that P^ does not end in v. 

V8 [Flip colours on Pi^] Recolour the edges on P^ by exchanging 0 and j. Down¬ 
shift from k, recolour f{vwk) - 0, and return to Step V2. ■ 

To see that this algorithm is correct, one needs to check that the recolourings 
in Steps V6 and V8 are legal. A careful analysis is given by Misra and Cries ||33l. 

For the running time, first note that Step V6 is repeated at most deg v times, 
so the algorithm eventually has to leave that step. The most time-consuming step 
is Step V7; a [0,/-path can be constructed in time 0{n) if for each vertext we 
maintain a table of incident edges indexed by colour. Thus the total running time 
of Algorithm V is 0{mn). 

Another example from this class of algorithms appears in the proof of Brooks’s 
theorem (see Chapter 2 and (SI), which relies on an algorithm that follows Algo¬ 
rithm G but attempts to re-colour the vertices of bichromatic components whenever 
a fresh colour is about to be introduced. 

Random changes 

There are many other graph colouring algorithms that fall under the umbrella of 
local transformations. Of particular interest are local search algorithms that re¬ 
colour individual vertices at random. This idea defines a random process on fhe sef 
of colourings called fhe Glauber or Metropolis dynamics, or fhe nafural Markov 
chain Monfe Carlo mefhod. The aim here is nol merely to find a colouring (since 
q > 4A, fhis would be easily done by Algorifhm G), buf fo find a colouring fhaf is 
uniformly disfribufed among all ^-colourings. 

Algorithm M {Metropolis) Given a graph G with maximum degree A and a q- 
colouring /q for q > 4A, this algorithm finds a uniform random q-colouring fj in 
polynomial time. 
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Ml [Outer loop] Set T - \qnln2n/{q-4A)']. Do Step M2 for t = 1,2,..., T, then 
terminate. 

M2 [Recolour a random vertex] Pick a vertex v e V and a colour c € [ 1,2,..., < 7 } 
uniformly at random. Set ft = ft-\. If c does not appear among v’s neighbours, 
then set ffv) = c. ■ 


An initial colouring /q can be provided in polynomial time because ^ > A + 1 
- for example, by Algorithm G. To see that the choice of initial colouring /o has 
no influence on the result fr, we consider two different initial colourings /o and 
and execute Algorithm M on both, using the same random choices for v and c in 
each step. 

Let dt = |(v: ffv) + //(v)}| be the number of disagreeing vertices after t 
executions of Step M2. Each step can change only a single vertex, so \dt - dt-\\ = 
1, 0, or -1. We have dt = dt-\ + 1 only if ft-i{v) - but /fv) + flG), 

so exactly one of the two processes rejects the colour change. In particular, v 
must have a (disagreeing) neighbour w with c - ft-\{w) //_i(w) or ft-i{w) 

fffw) = c. There are dt-i choices for w and therefore 2Adt-i choices for c and v. 
Similarly, we have dt = dt-\ - 1 only if ft-\{v) + ft-\{v) and c does not appear in 
v’s neighbourhood in either ft-\ or There are at least {q-2A)dt-\ such choices 
for c and v. 

Thus, the expected value of dt can be bounded as follows: 


,,(<?-2A)EW_i ] 2AEM_i] ^ - 4A 

L[dt\ < L[dt-\] +-= E[df-i] 1 - 


qn qn 

Iterating this argument and using do < n, we have 


qn 


, q-4A\^ I T{q-4A)\ , ^ , 

L{dT^ <n\l -) < n exp(---) < n exp(-In 2n) - j . 


qn 


qn 


By Markov’s inequality, and because dj is a non-negative integer, we conclude that 


Pr(/r = fr) = Pr(^^r = 0) > 1 - Pr(dr > 1) > 1 - E[r/r] > 5 ■ 


We content ourselves with this argument, which shows that the process is ‘suf¬ 
ficiently random’ in the sense of being memoryless. Informally, we can convince 
ourselves that fj is uniformly distributed because we can assume that /J in the 
above argument was sampled according to such a distribution. This intuition can 
be formalized using standard coupling arguments for Markov chains; our calcula¬ 
tions above show that the ‘mixing time’ of Algorithm M is 0{n log n). 

Algorithm M and its variants have been well studied, and the analysis can 
be much improved (see the survey of Frieze and Vigoda IfT^ l. Randomized lo¬ 
cal search has wide appeal across disciplines, including simulations in statistical 
physics and heuristic methods in combinatorial optimization. 
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6. Vector colouring 

We now turn to a variant of vertex-colouring that is particularly interesting from an 
algorithmic point of view. 


Vector chromatic number 


Let 5^“' = {x e ||x|| = 1). A vector q-colouring m d < n dimensions is a 
mapping x: V ^ 5^“' from the vertex-set to the set of r/-dimensional unit vectors 
for which neighbouring vectors are ‘far apart’, in the sense that their scalar product 
satisfies 

(x(v), x{w)} < -, for vw € E. 

q-1 

The smallest such number q is called the vector chromatic number x{G), which 
need not be an integer. For instance, the vertices of the 3-chromatic cycle graph Cj 
can be laid out on the unit circle in the form of a pentagram Then the angle 
between vectors corresponding to neighbouring vertices is ^tt, corresponding to 
the scalar product -1/(V5 - 1), so;y(C 5 ) < ^J5 < 3. 

Theorem 6.1 IfG has clique number a){G), then cjo{G) < x{G) < x(G). 


Proof. For the first inequality, let IT be a clique in G of size r = a>(G) and consider 
a vector ^-colouring x of G. Let y = liveiv x(v). Then 

0 < <y, y> < r • 1 -I- r(r - 1) • j, 

which implies that r < q. 

For the second inequality, place the vertices belonging to each colour class at 
the corners of a (q- l)-dimensional simplex. To be specific, let /: V {\,2,...,q} 
be an optimal ^-colouring and define x(v) = {x\,X 2 ,..., x„) by 


Xi = 


{{q-l)lqf''^, 

■ -{q{q - 1))“^^^ , 


0 , 


if i = /(v); 

if i + /(v) and i < q\ 

'\fi>q. 


Then we have 


(x(v), x(v)> 


q-1 q-1 

- --I- — - 

q q(q - 1) 


= 1 , 


and for v and w wifh f(v) + f{w) we have 


(x(v), x(w)) 



q- 1 

d 



q 

q-1 


) q(q - 1) 


1 

q- r 


1/2 
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Thus, ;c is a vector ^-colouring, soJ^(G) is at most q. ■ 

What makes vector colourings interesting from the algorithmic point of view is 
that they can be found in polynomial time, at least approximately, using algorithms 
based on semidefinite programming. The details behind those constructions lie far 
outside the scope of this chapter (see Gartner and Matousek Ifldll l. 

Theorem 6.2 Given a graph G wilhx{G) - q, a vector {q + e)-colouring ofG can 
be found in time polynomial in n and log(l/e). 

For a graph with a>(G) = xiG), Theorem 16. 1 1 shows that the vector chromatic 
number equals the chromatic number. In particular, it is an integer, and can be 
determined in polynomial time using Theorem l6.2l with £ < ^. This shows that the 
chromatic numbers of perfect graphs can be determined in polynomial time. The 
theory behind this result counts as one of the highlights of combinatorial optimiza¬ 
tion (see Grotschel, Lovasz and Schrijver ifThll l. 

How does the vector chromatic number behave for general graphs? For q = 2, 
the vectors have to point in exactly opposite directions. In particular, there can be 
only two vectors for each connected component, so vector 2-colouring is equivalent 
to 2-colouring. 

But already for q = 3, the situation becomes more interesting, since there exist 
vector 3-colourable graphs that are not 3-colourable. For instance, the Grotzsch 
graph, the smallest triangle-free graph with chromatic number 4, admits the vector 
3-colouring shown in Fig. [TO] as an embedding on the unit sphere. More com- 


V7 



V7 




Fig. 10: Left: the Grotzsch graph Middle and right: a vector 3-colouring 

plicated constructions (that we cannot visualize) show that there exist vector 3- 
colourable graphs with chromatic number at least (see ifT^ and |[22l ). 
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Randomized rounding 

Even though the gap between x and x can be large for graphs in general, vector 
colouring turns out to be a useful starting point for (standard) colouring. The next 
algorithm, due to Karger, Motwani and Sudan ll22Tl . translates a vector colouring 
into a (standard) vertex-colouring using random hyperplanes. 

Algorithm R {Randomized rounding of vector colouring) Given a 3-chromatic 
graph G with maximum degree A, this algorithm finds a q-colouring in polynomial 
time, where the expected size of q is E[g] = log n). 

R1 [Vector colour] Set e = 2 • 10~^ and compute a vector (3 -i- e)-colouring x of G 
using semidefinite programming. Let a > arccos(-l/(2 + e)) be the minimum 
angle in radians between neighbouring vertices. 

R2 [Round] Set 

^ = riOg;,/(;,_„)(2A)l 

and construct r random hyperplanes H\,H 2 ,..., //^ in R”. For each vertex v, 
let /(v) be the binary number brbr-\ ■ ■ - bi, where bi = 1 if and only if x{v) is 
on the positive side of the ith hyperplane Hi. 

R3 [Handle monochromatic edges recursively] Iterate over all edges to find the 
set of monochromatic edges M = {vw e E: f{v) - f{w) j. Recolour these 
vertices by running Algorithm R recursively on G[M], with fresh colours. ■ 

Figure [H] illustrates the behaviour of Algorithm R on the vector 3-colouring 
of the Grotzsch graph from Fig. [TOl Two hyperplanes separate the vertices into 
four parts. The resulting vertex-colouring with colours from [0,1)^ is shown to the 
right. In this example, the set M of monochromatic edges determined in Step M3 
contains only the single edge viovn, drawn bold in the figure. 



Fig. 11: Lefl and middle: fwo hyperplanes Righf: fhe corresponding colouring 


Algorifhm R algorifhm runs in polynomial lime, because Theorem 16.21 ensures 
lhal Step R1 can be performed in polynomial lime. 
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We proceed to analyze the size of the final colouring. Step R2 uses the colours 
{0,1,, 2'"“^), so the number of colours used in each Step R2 is 

2 '' < (2A)-l^ ^og(n/{n-a)) ^ , (5) 

what is more difficult is to bound the total number of recursive invocations. To this 
end, we need to understand how fast the instance size, determined by the size of M 
in Step R3, shrinks. 

Let e be an edge whose endpoints received the vector colours x and y. Elemen¬ 
tary geometrical considerations establish the following result. 

Theorem 6.3 Let x,y with angle tp (in radians). A random hyperplane in 
fails to separate x and y with probability 1 - (pjn. 

The angle between the vectors x and y is at most a. (To gain some intuition of 
this, if we ignore the error term e. Theorem 16.31 shows that x and y end up on the 
same side of a random hyperplane with probability 1 - o/tt < 1 - arccos(-^)/n' = 
1 - Inl^n = 5 -) The edge e is monochromatic if all r independent random hyper¬ 
planes fail to separate x and y in Step R2. Thus, 

Pr(e e M) < (1 - ajn)' < (nl(n - < 1/2A. 

By linearity of expectation, the expected size of M is 

E[|M|] = 2 Pr{e e M) < mjlA < \n. 

eeE 

Since each edge has two vertices, the expected number of vertices in the recursive 
instance G[M] is at most and in general, for i > 2, the expected number of 
vertices n,- in the fth instance satisfies n,- < In particular, n, < 1 after t - 

C?(log n) rounds, at which point the algorithm terminates. With the bound (l5]l on 
the number of colours used per round, we conclude that the total number of colours 
used is logn) in expectation. 

In terms of A, Algorithm R is much better than the bound of A -i- 1 guaranteed 
by Algorithm G. For an expression in terms of n, we are tempted to bound A by 
0{n), but that just shows that the number of colours is logn), which is 

worse than the 0{sjn) colours from Algorithm W. 

Instead, we employ a hybrid approach. Run Steps W1 and W2 as long as 
the maximum degree of the graph G is larger than some threshold d, and then 
colour the remaining graph using Algorithm R. The number of colours used by the 
combined algorithm is of the order of (In/d) + (Id)^'^^^ log n, which is minimized 
around d = with value 
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Variants of Algorithm R for general ^-colouring and with intricate rounding 
schemes have been investigated further (see Langberg’s survey 1271 ). The current 
best polynomial-time algorithm for colouring a 3-chromatic graph based on vector 
colouring uses colours, due to Chlamtac f9i 


7. Reductions 

The algorithms in this chapter are summarized in Tabled] 


Algorithm 

Time 

Problem 

B Bipartition 

0{n -1- m) 

2-colouring 

C Contraction 

0(1.619"+'") 

P(G,q) 

D Dynamic programming 

3"poly(n) 

XiG) 

G Greedy 

Oin + m) 

(A(G) -H l)-colouring 

I Inclusion-exclusion 

2" poly(n) 

XiG) 

L Lawler’s algorithm 

0(1.443") 

3-colouring 

M Metropolis dynamics 

poly(n) 

random q'-colouring (q > 4A) 

P Palette restriction 

1.5"poly(n) 

3-colouring 

R Rounded vector colouring 

poly(n) 

0(AO ®^i logn)-colouringforyf(G) _ ^ 

V Vizing’s algorithm 

0(mn) 

edge (A(G) -H l)-colouring 

W Wigderson’s algorithm 

0(n + m) 

0(-\/n)-colouring for;^f(G) = 3 

X Exhaustive search 

q'"poly(n) 

P(G,q) 


Table 1: Algorithms discussed in this survey 

Not only do these algorithms achieve different running times and quality guar¬ 
antees, they also differ in which specific problem they consider. Let us now be 
more precise about the variants of the graph colouring problem: 

Decision Given a graph G and an integer q, decide whether q can be (^-coloured. 

Chromatic number Given a graph G, compute the chromatic number ;t^(G). 

Construction Given a graph G and an integer q, construct a ( 7 -colouring of G. 

Counting Given a graph G and an integer q, compute the number P{G, q) of q- 
colourings of G. 

Sampling Given a graph G and an integer q, construct a random ^-colouring of G. 

Chromatic polynomial Given a graph G, compute the chromatic polynomial - that 
is, the coefficients of the integer polynomial q i-> P{G, q). 
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Some of these problems are related by using fairly straightforward reductions. 
For example, the decision problem is easily solved using the chromatic number 
by comparing q with ;^(G); conversely, xiG) can be determined by solving the 
decision problem for ^ = 1,2,..., n. It is also clear that if we can construct a q- 
colouring, then we can decide that one exists. What is perhaps less clear is the other 
direction. This is seen by a self-reduction that follows the contraction algorithm. 
Algorithm C. 

Reduction C {Constructing a colouring using a decision algorithm). Suppose that 
we have an algorithm that decides whether a given graph G can be ^-coloured. If 
G = K„ and n < q, give each vertex its own colour and terminate. Otherwise, select 
two non-adjacent vertices v and w in G. If G -i- vw cannot be ^-coloured, then every 
^-colouring / of G must have /(v) = /(w). Thus we can identify v and w and 
recursively find a (^-colouring for GI vw. Otherwise, there exists a (^-colouring of G 
with /(v) 4^ f{w), so we recursively find a colouring for G -i- vw. ■ 

Some of our algorithms work only for a specific fixed q, such as Algorithm B 
for 2-colourability or Algorithm L for 3-colourability. Clearly, they both reduce 
to the decision problem where q is part of the input. But what about the other 
direction? The answer turns out to depend strongly on q: the decision problem 
reduces to 3-colorability, but not to 2-colorability. 

Reduction L {q-colouring using 3-colouring). Given a graph G = {V, E) and an 
integer q, this reduction constructs a graph H that is 3-colourable with colours 
{0,1,2) if and only if G is (^-colourable with colours {1,2,...,^}. 

First, to fix some colour names, fhe graph H contains a triangle with the vertices 
0,1,2. We assume that vertex i has colour i, for / = 0,1,2. 

For each vertex v € V, the graph H contains 2q vertices vi,V 2 ,...,v^ and 
Vj, V 2 ,..., v^. Our intuition is that the v,s act as indicators for a colour in G in the 
following sense: if v; has colour 1 in // then v has colour / in G. The vertices are 
arranged as in Fig. [T^ al: the right-most vertex is 1 or 2, depending on the parity 
of q. The vertices vi, V 2 ,..., are all adjacent to 2, and so must be coloured 0 or 



2 - q (mod 2) 



1 


(a) 


(b) 


Fig. 12: 
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1. Moreover, at least one of them must be coloured 1, since otherwise, the colours 
for Vj, V 2 ,..., are forced to alternate as 1,2,1 ,..., conflicting with the colour of 
the right-most vertex. 

Now consider an edge vw in G. Let vi, V 2 ,..., and wi,>V 2 ,..., be the 
corresponding ‘indicator’ vertices in H. For each colour i = the vertices 

Vi and Wi are connected by a ‘fresh’ triangle as shown in Fig. fT^ bl. This ensures 
that V, and w, cannot both be 1. In other words, v and w cannot have received the 
same colour. ■ 

The above reduction, essentially due to Lovasz |[30l . can easily be extended to 
a larger, fixed q > 3, because G is ^-colourable if and only if G with an added 
‘apex’ vertex adjacent to all other vertices is {q + l)-colourable. For instance, 4- 
colourability is not easier than 3-colourability for general graphs. 

Thus, all (^-colouring problems for ^ > 3 are (in some sense) equally difficult. 
This is consistent with the fact that the case q = 2 admits a very fast algorithm 
(Algorithm B), whereas none of the others does. 

Many constructions have been published that show the computational difficulty 
of colouring for restricted classes of graphs. We will sketch an interesting example 
due to Stockmeyer If35l : the restriction of the case ^ = 3 to planar graphs. Consider 
the subgraph in Fig. HAl a). called a planarity gadget. One can check that this 



S 

(a) 



w 


V 


(b) 



w 


Fig. 13: A planarity gadget 

subgraph has the property that every 3-colouring / satisfies /(E) = /(W) and 
/(N) = /(S). Moreover, every partial assignment / to {N, S,E,W) that satisfies 
/(E) = /(W) and /(N) = /(S) can be extended to a 3-colouring of the entire 
subgraph. 

The gadget is used to transform a given (non-planar) graph G as follows. Draw 
G in the plane and for each edge vw replace each edge intersection by the planarity 
gadget. The outer vertices of neighbouring gadgets are identified, and v is iden¬ 
tified with W in its neighbouring gadget (see Eig. HAl bl). The resulting graph is 
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planar, and it can be checked that it is 3-chromatic if and only if G is 3-chromatic. 
Thus, the restriction to planar instances does not make 3-colourability computa¬ 
tionally easier. Unlike the case for non-planar graphs, this construction cannot be 
generalized to larger q > 3, since the decision problem for planar graphs and every 
q> 4 has answer ‘yes’ because of the four-colour theorem. 

Computational complexity 

The field of computational complexity relates algorithmic problems from various 
domains to one another in order to establish a notion of computational difficulty. 
The chromatic number problem was one of the first to be analysed in this fashion. 
The following reduction, essentially from the seminal paper of Karp ||23l, shows 
that computing the chromatic number is ‘hard for the complexity class NP’ by 
reducing from the NP-hard satisfiability problem for Boolean formulas on conjunc¬ 
tive normal form (CNF). This implies that ah other problem in the class NP reduce 
to the chromatic number. 

The input to CNF-Satisfiability is a Boolean formula consisting of s clauses 
Cl, C 2 ,..., Cs- Each clause Cj consists of a disjunction Cj - {lj\ V Zj 2 V • • • V IjC) of 
literals. Every literal is a variable xi,X 2 , ■ ■ ■ ,Xr or its negation xi,X 2 , ■ ■ ■ ,Xr. The 
problem is to hnd an assignment of the variables to ‘true’ and ‘false’ that makes all 
clauses true. 

Reduction K {Satisfiability using chromatic number). Given an instance Ci,C 2 , 

.. .,Cs of CNE-Satishability over the variables x\,X 2 ,..., Xr, this reduction con¬ 
structs a graph G on 3r - 1 - 5 - 1 - 1 vertices such that G can be coloured with r + \ 
colours if and only the instance is satisfiable. 

The graph G contains a complete subgraph on r - 1 - 1 vertices {0,1,..., r}. In 
any colouring, these vertices receive different colours, say /(/) = i. The intuition 
is that the colour 0 represents ‘false’, while the other colours represent ‘true’. Eor 
each variable Xj {I < i < r) the graph contains two adjacent ‘literal’ vertices v, 
and Vi, both adjacent to all ‘true colour’ vertices {1,2,..., r) except i. Thus, one 
of the two vertices v,-, v, must be assigned the ‘true’ colour i, and the other must be 
coloured 0. The construction is completed with ‘clause’ vertices wj, one for each 
clause Cj {I < j < s). Eet x,-,, ..., X 4 be the variables appearing (positively or 

negatively) in Cj. Then wj is adjacent to {0,1,..., r) \ {ii,i 2 , ■ ■ ■, 4)- This ensures 
that only the ‘true’ colours {/i, 4,..., 4) are available at wj. Eurthermore, if x, 
appears positive in Cj, then wj is adjacent to v,; if x,- appears negated in Cj, then 
Wj is adjacent to v,-. Eigure [14] shows the reduction for a small instance consisting 
of just the clause C\ = (xi V X 2 V X 3 ) and a valid colouring corresponding to the 
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assignment xi - xj - true, X 2 - false; the edges of the clique on {0,1,2,3) are not 
shown. Thus, the only colours available to wy are those chosen by its literals. ■ 



Fig. 14: A 4-colouring instance corresponding to C\ = {x\ V X2 V xj) 


Edge-colouring 

A mapping f:E {1,2,... ,^) is an edge-colouring of G if and only if it is a 
vertex-colouring of the line graph L{G) of G. In particular, every vertex-colouring 
algorithm can be used as an edge-colouring algorithm by running it on L{G). For 
instance. Algorithm I computes the chromatic index in time 2'”poly(n), which 
is the fastest currently known algorithm. Similarly, Algorithm G finds an edge¬ 
colouring with (2A - 1) colours, but this is worse than Algorithm V. In fact, since 
A < x'(G) < A -I- 1, Algorithm V determines the chromatic index within an addi¬ 
tive error of 1. However, deciding which of the two candidate values for;y'(G) is 
correct is an NP-hard problem, as shown by Holyer llT9l for;y'(G) = 3 and Leven 
and Gain 1291 for x'{G) > 3. 

Approximating the chromatic number 

Algorithm V shows that the chromatic index can be very well approximated. In 
contrast, approximating the chromatic number is much harder. In particular, it is 
NP-hard to 4-colour a 3-chromatic graph (see llTTl l. This rules out an approximate 
vertex-colouring algorithm with a performance guarantee as good as Algorithm V, 
but is far from explaining why the considerable machinery behind, say. Algorithm 
R results only in a colouring of size for 3-chromatic graphs. The best currently 
known exponent is c = 0.204 (see |[24l ). 

For sufficiently large fixed q, if is NP-hard fo find an exp(Q((7^^^))-colouring for 
a (^-colourable graph. If q is nof fixed, even sfronger hardness resulfs are known. 
We saw in Section 0fhaf fhe polynomial-fime compufable function ;y(G) is a lower 
bound on;if(G), even fhough fhe gap can sometimes be large, say ;y(G) > n^'^^''x{G) 
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for some graphs. Can we guarantee a corresponding upper bound for xl If not, 
maybe there is some other polynomial-time computable function g so that we can 
guarantee, for example, g{G) < xiO) < The answer turns out to be ‘no’ 

under standard complexity-theoretic assumptions: For every £ > 0, it is NP-hard 
to approximate;y(G) within a factor as shown by Zuckerman lldTI . 

Counting 

The problem of counting the (^-colourings is solved by evaluating P{G, q). Con¬ 
versely, because the chromatic polynomial has degree n, it can be interpolated using 
Lagrangian interpolation from the values of the counting problem at ^ = 0,1,..., n. 
Moreover, note that x^G) > q if and only if P{G, q) > 0, so it is NP-hard to count 
the number of (^-colourings simply because the decision problem is known to be 
hard. In fact, the counting problem is hard for Valiant’s counting class #P. 

On the other hand, an important result in counting complexity |[2T1l relates the 
estimation of the size of a finite set to the problem of uniformly sampling from it. 
In particular, a uniform sampler such as Algorithm M serves as a ‘fully polynomial 
randomized approximation scheme’ (FPRAS) for the number of colours. Thus, 
provided that q > 4A, Algorithm M can be used to compute a value g{G) for which 
(1 - e)g{G) < P{G, ^) < (1 -I- e)g{G) with high probability in time polynomial in n 
and 1 je for any £ > 0. Much better bounds on q are known (see the survey of Frieze 
and Vigoda |[T3]| '). Without some bound on q, such an FPRAS is unlikely to exist 
because, with s - 5 , it would constitute a randomized algorithm for the decision 
problem and would therefore imply that all of NP can be solved in randomized 
polynomial time. 


Conclusion 

Together, the algorithms and reductions presented in this survey give a picture of 
the computational aspects of graph colouring. For instance, 2-colouring admits 
a polynomial time algorithm, while 3-colouring does not. In the planar case, 4- 
colouring is trivial, but 3-colouring is not. An almost optimal edge-colouring can 
be found in polynomial time, but vertex-colouring is very difficult to approximate. 
If q is sufficiently large compared to A(G) then the set of colourings can be sampled 
and approximately counted, but not counted exactly. Finally, even the computation¬ 
ally hard colouring problems admit techniques that are much better than our initial 
Algorithm X. 

None of these insights is obvious from the definition of graph colouring, so the 
algorithmic perspective on chromatic graph theory has proved to be a fertile source 
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of questions with interesting answers. 
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